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SUMMARY 


A FORTRAN computer code for predicting the off-desigh performance of 
multistage axial-flow compressors is presented. The code, which was 
developed at the NASA Lewis Research Center, uses a meanline stage-stacking 
method. Stage and cumulative compressor performance is calculated from 
representative ineanlifiS velocity diagrams located at rotor inlet and outlet 
mean line radii. 

Numerous options are available within the code: 

(1) Nondimensional stage characteristics may be input directly or 
calculated froifl stage design performance input. 

{' 1 ) Stage characteristics may be modified for off-design speed and blade 
reset. 

(3) Rotor design deviation angle may be modified for off-design flow, 
speed, and blade setting angle. 

(4) Units of input and output may be SI or U.S. customary. 

Many of the code's options use correlations that are normally obtained from 
experimental data. These empirical correlations permit modeling the trends 
in stage and overall performance by a simple one-dimensional stage-stacking 
technique. However, the correlations may only be accurately applied to 
predict the performance of compressors similar to those compressors used in 
deriving the empirical correlations. The code is described in sufficient 
detail so that users may modify the correlations to suit their needs. 

Example calculations for a two-stage fan without blade reset and for 
three single stages with inlet-guide-varie reset agree well with experimental 
data. For off-design compressor performance prediction, the main features 
of the stage-stacking method are (1) simplicity, (2) fast convergence, and 
(3) the ability to directly incorporate correlations from experimental data 
to model real flow conditions. 


INTRODUCTION 

The axial-flow compressor is widely used in aircraft engines. In 
addition to its inherent advantage of high mass flow per frontal area, it 
can give very good aerodynamic performance. However, good aerodynamic 
performance over an acceptable range of operating conditions is not easily 
attained. A successful design and development process for multistage axial- 
flow compressors requires that numerous criteria be satisfactorily met. 

These criteria include (1) design optimization based on design and 
off-design performance considerations (ref, 1), (2) prediction of part-speed 
performance and assurance of part-speed stall margin, and (3) determination 
of required starting bleed and the amount of variability of the inlet-guide- 
vane (I6V) and stator-blade rows to match the stages. 

Both experimental and analytical programs can be used in the development 
process for multistage axial-flow compressors. Since compressor 
experimental development in test facilities is expensive and time consuming, 
any insight into the onset and location of troublesome flow regimes that can 
reduce the amount of teUing is very valuable. One natural information 
source is experimental data from similar compressor stagesi But such data 
in sufficient detail are rarely available. New compressors are usually 
extrapolations from the data of their predecessors. 

Analytical methods that contain good flow modeling are an alternative 
way of gaining the insight needed for compressor development. There are, of 
course, several levels of sophistication for analytical programs; but in 
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general, only the level of sophistication required to evaluate the relevant 
flow phenomenon is desired in order to minimize complexity and to give high 
computational efficiency. Compared with other more sophisticated two- and 
three-dimensional models for compressor flow, the stage-stacking method is 
very simple. The simplicity of a one-dimensiohal compressible flow model 
enables the stage-stacking method to have excellent convergence properties 
and short computer run time (ref. 2). The simplicity of the model results 
in manageable computer codes that ease the incorporation of correlations 
directly linked to experimental test data to directly model real flow 
phenomena. 

The stage-stacking computer code discussed in this report was developed 
and used at the NASA Lewis Center during the past several years. It has 
been routinely used to generate performance maps for compressors evaluated 
experimentally at Lewis. Correlations from experimental data to model real 
flow phenomena were added so that the code's performance predictions agreed 
with the measured performance. The code is an extension of the stage- 
stacking method of reference 3. The present code either accepts 
nondifflensional stage characteristics as input or will calculate these 
characteristics from aerodynamic input available from compressor design 
codes such as reference 4. 

The code is described in sufficient detail herein to permit a user to 
modify the correlations from experimental data within the code. It is 
anticipated that at times revised correlations may better suit the 
particular needs of the user. The code itself is written in FORTRAN, and 
U.S. customary units are used in the coded correlations and calculations. 


CALCULATION PROCEDURE 

The calculation procedure is discussed in three parts. First, the 
stage-stacking method is described. Second, optional calculations available 
to the computer code user are described. Third, the computer code STGSTK is 
outlined. _ 
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Figure Z - Typical meanllne Inlet and outlet velocity diagrams for a rotor. 


Stage-Stacking Method 

To describe the stage-stacking method, first the flow assumptions are 
discussed and then the stage characteristics and the stacking procedure. 

Flow assumptions . - One-dimensional compressible flow is assumed. Flow 
continuity can therefore be expressed as 

W = pAV^ (1) 

where A is the annulus area. All symbols are defined in appendix A. 

Flow continuity is satisfied in the axial velocity V£ calculation at 
the rotor inlet and outlet axial locations (fig. 1) of each stage. Thus, 
for a given flow and speed and stage inlet flow conditions of total pressure 
ana temperature and absolute flow angle B2M> ® meanline velocity diagram 
(fig. 2) can be obtained at the rotor inlet. And, by assuming that the 
stage overall pressure ratio and adiabatic efficiency apply at the rotor 
exit, a meanline velocity diagram can be obtained at the rotor exit. If 
rotor exit total pressure and temperature are then assumed to apply at the 
inlet of the next rotor, meanline velocity diagrams can be obtained at every 
rotor inlet and exit from the overall stage performance parameters of 
pressure ratio Pr and adiabatic efficiency nad* These rotor inlet and 
outlet meanline velocity diagrams obtained from overall stage performance 
parameters are assumed to represent the stage and are referred to as 
representative meanline velocity diagrams within this report. Alterations 
to these representative meanline velocity diagrams are assumee to alter the 
associated stage performance parameters. Specific calculations are used to 
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vary the representative meanline velocity diagrams and to predict changes in 
the associated stage performance parameters Pr and na(j. 

Stage characteristics . - The stage performance characteristics consist 

of three nondimensional quantities - adiabatic efficiency nad» pressure 
coefficient and flow coefficient <p. They are calculated from 


S^ad (T3 - T^) 


'ad '‘3 
,.2 


Jy-D/’ 


The stage characteristics are usually presented as adiabatic efficiency 
versus flow coefficient nad(‘i’) and pressure coefficient versus flow 
coefficient >!'(?). Figure 3 shows typical stage characteristics for 
which the stage design point (reference condition) is at peak adiabatic 


Desigaicfitence) point 


Slige flow coefficient, (p 

Figure 3. - Twicit nondimensional stage characteristics. 



epiciency. At some low flow coefficient (/c the stage will stall, and 
at some high flow coefficient the stage will choke. Ideally the 
stage characteristics are independent of both compressor size and speed. 

One option of the computer code STGSTK is either to input the stage 
characteristics nad(<p) and *(<p) or to input only the overall stage 
pertormance Pr and at a design or reference point and have the 

computer code calculate the stage characteristics. If the code option is 
used to input only the overall stage performance Pr and at a 

design or reference point, the calculated stage characteristics are obtained 
from representative meanline velocity diagrams at the rotor inlet and 
outlet. These velocity diagrams are used to represent the overall stage 
flow, pressure ratio, and adiabatic efficiency and not just blade element 
performance at the midspan point. 

Stackin g procedure . - Once the appropriate stage characteristics are 
obtained, the stage-stacking procedure involves a straightforward 
calculation process. The overall compressor inlet flow conditions must be 
known, and they are used as the overall inlet flow condition for the first 
stage. Then, for various selected compressor speeds and flows, the 
calculation process is repeated for each stage as follows: 

( 1 ) Calculate the representative meanline velocity diagram at the rotor 
inlet and then the stage flow coefficient 9. 

(2) From the stage characteristics and obtain 

the stage overall adiabatic efficiency nad pressure coefficient. 

( 3 ) Calculate the representative meanlme velocity diagram at the rotor 
outlet and then the overall total pressure P and temperature T at the 
stage outlet, and use these values for the next stage inlet P and T. 

This process is repeated for each stage; and the cumulative compressor 
performance, which consists of compressor overall adiabatic efficiency 
temperature ratio, and pressure ratio, is calculated. The end result \s a 
compressor overall performance map of adiabatic efficiency and pressure 
ratio versus flow for various speeds. 


Optional Calculations 

Optional calculations performed within the computer code STGSTK are 
executed for three major reasons; code input selection, blade reset 
conditions, and stage characteristic adjustments for real flow effects. 
Concerning optional calculations for code input selection, the primary 
9Ption is whether or not stage characteristics nad(<p) and 4>(<p) are 
input. This input option has been discussed in the section Stage 
characteristics . If stage characteristics and 41(9) are not 

input, the code will calculate them based on the input of overall stage 
performance Pr and nad ^ design or reference point. 

For blade reset conditions, optional calculations are performed within 
the STGSTK code to alter the stage characteristic 9(9). An example of 
calculated changes in the stage characteristic 9(9) because of reset of 
upstream stator vanes is shown in figure 4 . Details of how the code STGSTK 
alters the stage characteristics for blade reset aYn of a blade with setting 
angle yq are given in appendix B. Basically the blaoe reset AYn alters 
the flow angles of the meanline representative velocity diagrams associated 
with the stage. New rep»"esentati ve velocity diagrams are calculated that 
determine the new stage characteristic 9(9) for the stage. 
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Figure 4. - Example of calculated changes of stage characteristic for a reset 
of upstream stator vanes. 


The remaining optional calculations adjust the stage characteristics for 
real flow effects that are not directly modeled by one-dimensional 
compressible flow. Three of these optional real flow adjustments are 
available. The first real flow adjustment alters the stage characteristic 
41(9) for part-speed conditions. This option is especially applicable 
for stages that have high inlet relative Mach numbers. Figure 5 shows the 
nature of the optional stage characteristic 4)(<p) adjustment for part 
speed. At the part-speed condition, as compared wi h design speed, the 
pressure coefficient 'll drops by an amount Avf and the range of the 
flow coefficient <p is expanded. 

The second real flow adjustment option alters the stage characteristic 
nad('p) part-speed conditions. This adjustment option for nad(<p) 
at part speed consists of two parts: (1) the expansion of the range of the 

flow coefficient <p, which is identical to the adjustment used for the 
stage characteristic 'l»(«p) discussed in the previous paragraph, and 
(2) a change in the level of nad at part speed, which is controlled by 
values of the input variable ETARAT, which is discussed later in the section 

Input Data. 

The third real flow adjustment option involves alterations of rotor 
deviation angle 6 r under any combination of three conditions: off-design 

flow coefficient «p, rotative speed N, and blade setting angle yq* ^ 
change in the, rotor deviation angle 6 r changes the rotor outlet relative 
flow angle 63M* This in turn alters the rotor outlet representative 
meanline velocity diagram and changes the stage characteristic 


Computer Code Outline 

The computer code STGSTK consists of a main routine and eight 
subroutines* Figure 6 gives a line representation of the subroutine calls. 
This section briefly summarizes the calculations within each subroutine. 
Details of the calculations are given in appendix B. 

The main computer routine is entitled MAIN* MAIN is the central control 
routine and it calls the major subroutines. The major subroutines do three 
things: (1) process the input arid output data, (2) calculate parameters for 

the compressor design (reference) point, and ( 3 ) oerform the optional 
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Figure 6. * Line represenatloil of jubroutiiifi calls for stage-stacking program. 


calculations discussed previously. A list of the subroutines and their 
primary purposes follows. 

CSINPT - read and write the input data 

CSPREF - calculate design (reference) parameters 

CSETA - obtain optional stage characteristic nad(«p) at design speed 

CSPSI - obtain optional stage characteristic «p(<p) at design speed 

CSPSD - perform option to alter pressure coefficient at off-design 
speeds 

CSPAN - perform option to alter stage characteristic 4 i(<p} because of 
blade reset 

CSOUPT - calculate and write the output data 

CPF - calculate specific heat Cp ana its ratio y 

Other optional calculations are located within the STGSTK code as 
follows. MAIN contains an option to alter the flow coefficient 9 at 
off-design speeds, CSPSI, CSPSD, and CSPAN contain an option to alter the 
rotor deviation angle 6 r for off-design flow, speed, and blade setting 
angle, res'*- ctively. 


COMPUTER CODE USER INFORMATION 

This section provides information for someone who wants to use the 
STGSTK code. The input data are described and an example input data set is 
given. The output data computed by STGSTK are also described. For a guide 
to the optional calculations that may be selected, the user may refer to the 
section Optional Calculations presented previously. 


Input Data 

All input data needed to use the stage-stacking program are described in 
this section. The input data are described in the same sequence as they are 
used by the program. Except for the title card and specific-heat polynomial 
coefficients CPCO, all the input format is for floating-point numbers in 
fields of 10. Figure 7 shows the location of input. Except for CPCO the 
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Fiqure 7. - Input variable locations on cards foi* stage-stacking program. 


units of the input data can be selected as either all SI or all U.S. 
customary. Program subroutine CSINPT reads and prints the data. 

The following is a list of the input data as they are read by subroutine 
CSINPT. For each input variable name listed, its format, description, and 
units are included. 


TITLE, 18A4 

title card on which any alphanumeric data can be used; one 
card needed 

STAGEN, 

Flo 

number of stages 

SPEEDN, 

FIO 

number of speed lines 

CHAPTS, 

FIO 

number of points used to describe stage characteristic 

PO, FIO 


inlet total pressure, N/cm^ (psi) 

TO, Flo 


inlet total temperature, K (“R) 

WTMOLE, 

FIO 

molecular weight 

DESRPM, 

FIO 

design rotative speed, rpm 

DESFLO, 

FIO 

design flow, kg/sec (Ib/sec) 

SPDPSI, 

Flo 

alters 4> value for off-design speed when equal to 1.0 

SPDPHI, 

Flo 

alters 9 value for off-design speed when equal to 1.0 

DRDEVG, 

Flo 

alters rotor deviation angle for blade reset when equal to 
1.0 


8 


II iliiPiii I 


DRUEVN, Flo 

DRDEVP, FIO 

UNITS, Flo 

CPCO, E20.8 

RT2, FIO 
RH2, FIO 
RT3, FIO 
RH3, Flo 
BET2M, FIO 
CB2M, FIO 

CB2MR, FIO 

CB3MR, FIO 

RK 2 M, Flo 
RSOLM, Flo 
SK2M, FIO 
PR, Flo 
ETAINP, FIO 
PCTSPD, Flo 

ETARAT, Flo 

BLEED, Flo 
PHIDES, FIQ 


alters rotor deviation angle for off-design speed when 
equal to 1.0 

alters rotor deviation angle for off-design when equal 
to 1.0 

used to specify units of input; use 1.0 for 31, 0.0 for 
U.S. customary 

specific-heat Cp polynomial coefficients in U.S. 
customary units (Btu ’R-l/lb ... Btu "R^^/lb) 

rotor inlet tip radius, cm (itu) 

rotor inlet hub radius, cm (in.) 

rotor outlet tip radius, cm (in.) 

rotor outlet hub radius, cm (in.) 

rotor inlet absolute flow angle at meanline radius, deg 

change in rotor inlet absolute flow angle atmear.line 
radius, deg 

change in rotor inlet relative flow angle at meanline 
radius, deg 

change in rotor outlet relative flow angle at meanline 
radius, deg 

rotor inlet blade metal angle at meanline radius, deg 

rotor blade row solidity at meanline radius 

stator inlet blade metal angle at meanline radius, deg 

design stage pressure ratio used to calculate PSIDES 

design stage adiabatic efficiency used to calculate ETADES 

value of rotative speed expressed as a decimal fraction of 
design speed; design speed value or 1.0 must be the first 
value for this input variable 

ratio of adiabatic efficiency at design speed to adiabatic 
efficiency at speed corresponding to PCTSPD; 1.0 is 
normally the first value for this input variable 

bleed flow for a particular stage and speed corresponding 
to PCTSPD, kg/sec (Ib/sec) 

stage flow coefficient at design speed 
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PSIDES, Flo 

ETADES, Flo 

SPEEDF, FIO 
FLOWIN, FIO 

DFLOW, FIO 

FLOWFI, FIO 


stage pressure coefficient at design speed; when input PR 
is not zero, PSIDES must be zero 

stage adiabatic efficiency at design speed; when input 
ETAINP is not zero, ETADES must be zero 

decimal fraction of design speed for a particular speed line 

value of lowest flow for speed line designated by 
SPEEDF, kg/sec (Ib/sec) 

change in flow for speed line designated by SPEEDF, 
kg/sec (Ib/sec) 

value of highest flow for speed line designated by 
SPEEDF, kg/sec (Ib/sec) 


Example Input Data Set 

The example program data set listed in this section is for the NASA 
Lewis two-stage fan having the low-aspect-ratio, first-stage rotor blading 
of reference 5. The fan design information is used for the input data, bi 
units are used except for the Cp polynomial coefficients. 


NASA TWO 
2. 

1 . 

0, 
0 . 

25.A20 
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24.628 
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. 0 
.0 
.0 
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.0 
1 . 
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.8 
. 7 
. 5 
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5. 

1 . 

23746571E 00 
13991136E-09 
23.960 
13.604 
23.566 
14.696 
.0 
.0 
.0 
. 0 

55.46 

1.57 

36.15 

1.5087 

.8690 

.9 

1.017 


,35 


.42 


32. 
27 . 
21 . 
17 . 
11 . 


LOW R1 AR, TP-1493 
8. 10.135 

1 . 1 . 

0.219619996-04 
-0. 78056154E-13 


288.17 

1 . 

-0 

0 


51 INPUT 
28.97 
1 . 

.87791479E-07 
. 15042604E-16 


.8 

.7 

.5 


1.029 

1.017 

1.023 


.38 

.42 

.43 

.44 

.44 

.45 

.46 

.48 

.2 

35. 



.5 

32. 



.5 

28. 



.5 

24. 



.5 

18. 
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Output Data 

All output data generated by the stage-stacking program are described in 
this section. The output data consist of two types: 

(1) Intermediate output for changes in the stage characteristics because 
of blade reset and off-design flow conditions 

(2) Final output for the predicted compressor performance along the 
selected speed lines 

The program's MAIN routine writes the intermediate output. Subroutine 
CSOUPT writes the final output. The units of the output data are either all 
SI or all U.S. customary depending on the value of the input parameter 
UNITS. For SI output, UNITS equals 1.0; for U»S. customary output, UNITS 
equals 0.0. 

The following is a list of the output data in the same order as they are 
written by the program. For each output variable lisced, its r.scription 
and units ahe included. .1 


PHIREF 

PSIREF 

DPHIA 

DPSIA 

F LOCAL 

BET2M 

BET3MR 

RINCM 

RDFM 

SINCM 

DPSIS 

PHI 


Stage flow coefficient at design flow, speedy and blade 
setting angles 

stage pressure coefficient at design flow, speed, and blade 
setting angles 

change in flow coefficient PHIREF because of rotor- or 
stator-blade reset 

change in pressure coefficient PSIREF because of rotor- or 
stator-blade reset 

calculated flow for a stage at design speed if blade reset 
is specified, kg/sec (Ib/sec) 

rotor inlet absolute flow angle at design speed and flow 
and specified blade reset, deg 

rotor outlet relative flow angle at design speed and flow 
and specified blade reset, deg 

rotor incidence angle at design speed and flow and 
specified blade reset, deg 

rotor diffusion factor at design speed and flow and 
specified blade reset, deg 

stator incidence angle at design speed and flow and 
specified blade reset, deg 

change in pressure coefficient PSIREF because of off-design 
speed effects 

calculated stage flow coefficient for selected speeds with 
blade reset effects included 
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PSI calculated stage pressure coefficient for selected speeds 

with blade reset effects included 

eta calculated stage adiabatic efficiency for selected speeds 

TAU stage temperature rise coefficient 

PP stage pressure ratio 

C-ETA compressor cumulative adiabatic efficiency 


EXPERIMENTAL DATA COMPARED WITH CODE PREDICTIONS 

The example calculations by the stage-stacking code discussed in this 
section are of two types: (1) calculations to predict overall compressor 

performance, and (2) calculations to predict flow reduction by inlet vane 
reset. The example calculations illustrate the procedure used for the 
stage-stacking code. The code's predictions are compared with experimental 
test data. 


Two-Stage Fan Performance 

These calculations were performed for the two-stage fan of reference 5, 
which has a design pressure ratio of 2.4 and an inlet tip speed of 429 
m/sec. The previous section Example Input Data Set lists the input data 
set. These input data select the following program calculation options: 

(1) stage characteristics nad(<c) and <|«((p) will be calculated from 
a reference design-point stage pressure Pr and adiabatic efficiency 
nad» (2) 4«(<p) will be altered for off-design speed, and (3) rotor 
deviation will be adjusted for off-design speed and flow. Appendix C lists 
a.11 the program output data. Stage performance and cumulative compressor 
predicted performance are listed for various flows at five selected speeds. 

In figure 8 the overall fan performance calculated by the stage-stacking 
program is compared with the experimental measured performance reported in 
reference 5. At design speed for any given pressure ratio the measured flow 
is greater than the calculated flow. This occurs because the calculated 
performance at design speed was forced throuah the fan's design point, and 
the measured data indicated that the fan performed at a higher flow than its 
intended design flow. For the off-design part speeds the discrepancies 
between the calculated and experimental data are similar but less severe. 
This would tend to support the credibility of the program to predict overall 
compressor performance at part-speed, off-design flow. No judgment can be 
made from these fan data on the program's ability to account for vane reset 
since the fan was tested with no inlet guide vanes and with fixed stators. 


Single-Stage Performance with Variable Inlet Guide Vanes 

When an inlet guide vane (IGV) is reset to increase the absolute flow 
angle ep at the following rotor inlet, the corrected flow at the rotor 
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Figure 8. - Comparison of calculated performance by stage-stacking 
code to experimental data for NASA two-stage fan of reference 5. 


inlet will be reduced. Also, the stage pressure ratio Pr will decrease, 
and this may cause the corrected flow at the stage exit to decrease. 
Corrected-f low reduction at the stage exit will depend on the amount of IGV 
reset, the design speed, and the stage performance. 

Calculations are performed for each of the three single-stage 
compressors with variable inlet guide vanes for which experimental data were 
reported in references 6 to 8. For these stages, the design rotor inlet tip 
speeds are 347, 457, and 427 m/sec, and the measured stage pressure ratios 
at peak stage efficiency are 1.42, 1.72, and 1.52, respectively. For each 
stage a comparison of the calculated to the experimental overall performance 
data at various IGV setting angles is shown in figures 9 to 11. Figure 12 
shows the corrected-f low reduction ratio at the stage outlet versus IGV 
setting angle at peak stage efficiency for design speed and 80 percent of 
design speed. The open symbols are for measured data and the lines are for 
calculated output from the stage-stacking program. Figure 13, which was 
derived from fiqure 12, shows the IGV reset angle required to reduce the 
corrected flow (at the stage outlet) by 10 percent versus rotor inlet tip 
design speed for 100 and 80 percent of design speed. 

Agreement between the calculated and measured IGV reset for a given flow 
reduction and tip speed is better at design speed than at 80 percent of 
design speed. However, for both speeds the trends of the calculated and 
measur >d data are very similar. This indicates that as the rotor speeo goes 
up, IGV reset must be increased for a given flow reduction ratio W/Wj. 

This influence of rotor speed on the vane-reset flow reduction relationship 
was previously discussed in reference 9. 
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(a) IGV setting angle. 0°. (b) IGV setting angle. 10®. (c) IGV setting angle. 20®. 

Figure 9. ♦ Comparison of calculated to experimental performance data for the 347-m/ sec-design-speed single-stage com- 
pressor of reference 6 at various inlet guide vane setting angles. 
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Figure 10. - Comparison of calculated to experimental perfor trance data for the 45f-mfsec -design -speed single-stage com- 
pressor ot reference ? at various inlet guide vane setting angles. 
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(b) 80 Percent of design speed. 

figure IZ - Corrected flow reduction ratio it stage outlet, at peak effi- 
ciency, as a function of Inlet guide vane setting angle tor three 
single-stage compressors. 
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(a) 80 Percent of design speed. (b) Design speed. 


Figure 13. - Required inlet guide vane reset for a 10-percent 
corrected flow reduction at the stage outlet, at peak efficiency, 
as a function of rotor in'et tip speW. 


CONCLUDING REMARKS 

A computer code for predicting off-design performance for multistage 
axial-flow compressors has been discussed in this report. The meanline 
stage-stacking method used has the following properties: 

(1) It is a one-dimensional, compressible flow model with fast 
convergence. 

(2) Overall stage performance is represented by meanline velocity 
diagrams at the rotor inlet and outlet. 

(3) Options are included to calculate the stage characteristics and to 
adjust them for blade reset and real flow effects. 

(4) Experimental test cata can be applied directly to correlations of 
model real flow conditions. 

(b) Accurate off-design predictions can be made for a limited range of 
compressors. 

Example calculations compared with experimental data for the stage- 
stacking code reported herein give the following indications: 






( 1 ) The code's calculation options to alter stage characteristics and 
rotor deviati^angles for off-design conditions resulted in a Performance 
preS^ctlon f7a two!stage, 2.4-pressure-ratio fan that compared well with 

experimental data ^ characteristics and rotor 

deviation angles L^to blade reset resulted in a flow reduction prediction 
with inlet-guide-vane reset for three single-stage compressors that compare 
satisfactorily with experimental data trends. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, October 8, 1981 


I U lliWPI 




APPENDIX A 
SYMBOLS 


The following is a list of symbols defined as used in the text and 
figures of this report. Only basic SI units for these symbols are given. 
The user information section of the report defines the computer program 
input and output data with appropriate SI and U.S. customary units as 
required by the program. 


A 

Cp 

••• 

Df 

i 

N 

P 

Pr 

R 

T 

Tr 

t 

U 

V 


annulus area, m^ 

specific heat at constant pressure, J/kg K 

coefficients for Cp polynomial, JK~l/kg ... JK~5/kg 

diffusion factor 

incidence angle, deg 

rotative speed, rpm 

total pressure, N/m^ 

pressure ratio 

gas constant, d/kg K 

total temperature, K 

temperature ratio 

static temperature, K 

wheel speed, m/sec 

velocity, m/sec 


W 

e 

Y 

YQ 

AYO 

6 

"lad 

K 

P 


flow, kg/sec 
flow angle, deg 
ratio of specific heats 
blade setting angle 
blade reset 
deviation angle, deg 
adiabatic efficiency 
blade metal angle, deg 
static density, kg/m^ 


0 

blade solidity 

9 

flow coefficient 

9 

pressure coefficient 

Subscripts: 

c 

choke 

d 

design condition 

H 

highcsc value 

h 

hub 

i 

indicates <p, N, or 

L 

lowest value 

M 

mean line 

N 

speed 

R 

rotor 

S 

stator 

s 

stall 

T 

tangential 

t 

tip 

Z 

axial 

9 

flow coefficient 

2 

rotor inlet 

3 

rotor outlet 

Superscript: 

( )• 

relative to rotor 


subscript 
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CODE CALCULATION DETAILS 

Details of calculations within the "^‘'“^.ra'lSsecU^ for the 

Generally a routine's common to several routines and program 

Srthfp?o^ra™?“‘A"guidftr?^^rsroSe°ce is the calculation flow chart of 
figure 14. 


Main Routine 

The main computer “de routine is entitled^MMN.^and^lt^ca^ 

3futinrfairs'rrom"wrN ard brief descriptions of the calculations 


of the 
within 
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each subroutine called from MAIN will be used to discuss the overall program 
structure. Other major items performed within MAIN are also discussed. 

MAIN first calls subroutine CSINPT. The primary purpose of subroutine 
CSINPT is to read and write thd input data required to run the computer 
code. After calling CSINPT, MAIN calculates several parameters associated 
with the rotor inlet and outlet for each stage. Among these parameters arr. 
the rotor inlet and outlet raeanline radii, which are calculated from 



Symbols are defined in appendix A. At these meanline radii, representative 
meanline velocity diagrams will be calculated for each stage. 

MAIN next calls subroutine CSPREF, which calculates design (reference) 
velocity diagrams at the meanline radii for each stage. The design velocity 
diagrams are based on the design values of speed, flow, and blade setting 
angles. CSPREF also calculates the design values of (reference) flow 
coefficient <pd, pressure coefficient 4>ij, and adiabatic efficiency 
Had d ^0^ stage from the following general expressions: 

. . ^ (B3) 

^2T 


= 


'ad 


S'^ad ^ '3 ‘2 


(T3 - T^) 

Pf,(Y-l)/Y _ I 
— 


(B4) 


(B5) 


The values of Cp and y are a function of static temperature and are 
obtained from subroutine CPF, whenever needed, by all subroutines throughout 
the calculation procedure. 

If the stage characteristic nad(<p) design speed is not input, 

MAIN calls subroutine CSETA to obtain nad(«p) f®'" stage at design 
speed. If the stage characteristic »ii((p) at design speed is not input, 

MAIN calls subroutine CSPSI to obtain for the stage at design 

speed. Also, MAIN has an optional call to subroutine CSPSD, which alters 
the pressure coefficient values for off-design speeds. If a rotor or stator 
is reset, MAIN calls subroutine CSPAN, which alters the stage 
characteristic <i«((fi) affected by the reset. Within MAIN there is an 
optional calculation that alters the flow coefficient tp values for 
off-design speeds. Within subroutines CSPSI, CSPSD, and CSPAN there is an 
optional calculation that alters rotor deviation angle for off-design values 
of flow, speedy and blade setting angle, respectively. 

MAIN writes the intermediate output data, which consists primarily of 
the following: (1) design values for each stage 9, <i, and nad 
specified flow angles, (2) changes in design values of stage 9 and 9 
because of blade reset, and (3) tabulated stage characteristics of 9 ( 9 ) 
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and nad(‘p) at specified speeds. These stage characteristics are used 
to calculate the compressor performance map. 

The final subroutine called by MAIN is CSUUPT. Subroutine CSOUPT reads 
the selected speeds and flows at which compressor performance is desired. 

At the selected speed and flow conditions CSOUPT calculates and writes 
individual stage and cumulative compressor performance parameters. 


Subroutine CSINPT 

The primary purpose of subroutine CSINPT is to read and write the input 
data that the program requires. The main text of this report has a section 
entitled Input Data that contains the information that a user needs to 
prepare an input data set for this program. This input data section 
explains the input data format, purpose, and units. CSINPT is coded with an 
option to enable the units of the input to be either all 51 or all U.S. 
customary. The write statements within CSINPT do not contain units and are 
therefore applicable to both SI or U.S. customary input data. CSINPT writes 
the input in the same units in which the input was read. A portion of the 
final lines of coding of CSINPT converts input of SI units into U.S. 
customary units. This conversion of the units of SI input data into U.S. 
customary units permits the FORTRAN expressions for the program calculations 
to be formulated in terms of U.S. customary units. 

For input of design stage performance there is an option of either of 
the following inputs: (1) stage pressure ratio Pr and adiabatic 

efficiency nad> (2) stage characteristics, which consist of pressure 
coefficient versus flow coefficient ^{(f) and adiabatic efficiency 
versus flow coefficient nad(‘p)* When either of these two input options 
is used as input, the input parameters for the option not used are input as 
zeros. Program subroutines called after CSINPT calculate values for the 
parameters that were not input by option. 


Subroutine CSPREF 

At design speed and flow, subroutine CSPREF calculates (1) velocity 
diagrams at the meanline radii for each rotor inlet and outlet, and ( 2 ) 
selected performance parameters for each stage. Figure 2 shows the meanline 
velocity diagrams associated with a typical rotor. CSPREF performs a 
one-dimensional, compressible, inviscid flow calculation at each rotor inlet 
and outlet to obtain the meanline velocity diagrams for design input 
conditions. 

The sequence of calculations in CSPREF is as follows: 

(1) From input design W, N, and A^, calculate by iteration 

(with Cp and y functions of t^) the design values for V 22 M and 9 . 

(2) If design stage Pr and naq input, calculate design i|). 

(3) If stage characteristics vw) and na(j{ip) at design speed are 
input, obtain t, and Pr by linear interpolation. 

(4| Calculate e^M and Vj2M* 

( 5 ) For design N and W, calculate by iteration the design values for 
Vt 3 m and V 23 m» with Euler's equation solved for Vj 3 m as 


and C. 

j f ^ 

S3M, the rotor and stator incidence angles, 

and the rotor diffusion factor by using 


Op and Y functions of t3« 
(6) Calculate design ' 


MR “ 

®2M 

" ''2MR 

(B7) 

MS " 

^3M 

“ '^2MS 

(B8) 

1 _ 

1 - 

1 

* '^3 m''t 3 M " '^2m''t2M 

(B9) 

fR 

'^2M ^*^3M '^2M^'^Rm'^2M 



These calculations within CSPREF are repeated for each stage of the 
compressor for which input was read by CSINPT. 


Subroutine CPF 

This subroutine is used to obtain values of Cp and y as a function 
of static temperature t. CPF is called by the previously discussed 
subroutine CSPREF and also is called by subroutines CSPSI, CSPSD, CSPAN, and 
CSOUPT. Subroutine CPF calculates Cp from a fifth-degree polynomial of 
t expressed by 

Cp » C^ + C2t + C3t^ + C^t^ + Cgt^ + Cgt^ (BIO) 

where the polynomial coefficients Cj to C5 are input data read by 
CSINPT. The value of y is then calculated from 



CPF also calculates various other functions of y used by the calling 
subroutines for the flow calculations. 


Subroutine CSETA 

This subroutine is called by MAIN when values of stage characteristic 
nad(<p) at design speed are not usable input (i.e., the input value for 
^ad. ■'5 CSllIK obtains values of ngd for each stage at the 

various input 9 for the stage. The following procedure is used within 
subroutine CSETA: 

(1) A curve is generated for nad('f) as depicted in figure 15. This 
curve consists of two parabolas. The first parabola extenas from the 

minimum flow (stall) coefficent 95 to the design flow coefficient 
9(j. The second parabola extends from the design flow coefficient 
9d to the maximum flow (choke) coefficient 9^. 

(2) For each stage the generated curve for nad(<p) •''as the 
properties: 
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Figure 15. *• Properties ot curve fit for stage adiabatic efficiency as a function 
of stage flow coefficient generated by subroutine CSFTA. 



“'’ad 

H = 0 

dq) 

at 

9 = 9^ 

’’ad ■ '’adjd 

at 

9 = 9d 

”ad = ”ad,d 

at 

9 = 9s 

’’ad = '’ad.d 

at 

9 = 9c 


The peak nad is located at the design (reference) condition. 

(3) After the curve as described above is generated for stage 

nad(<p)» subroutine CSETA then calculates a nad value for every 
input <p for the stage. 

This procedure is repeated within CSETA for each stage of the 
compressor. The merit of this particular procedure depends on its ability 
to simulate performance for the type of compressor being studied. Another 
procedure may better meet the needs of the user. The isolation of this 
procedure within a single subroutine readily permits its identification for 
alteration. 


Subroutine CSPSI 

This subroutine is called by MAIN when values of stage characteristic 
ip(a) at design speed are not usable input (i.e., the input value for 
^ is 0.0). CSPSI obtains values of for each stage at the various 




) 


j 

] 

! 

j 

I 
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input <p for the stage. The following calculdtion procedure is usr>d 
within subroutine CSPSI: 

(1) Assume design values calculated in CSPREF for e?MH. e?wn, 

9d. and 03(vic|. 

(2) For an input 9, calculate V22M froitl 


'^Z2M = '^Z2Md 


fe) 


calculate the flow W corresponding to the input 9 from 


Wo = 


^2^z'^Z2M 


value for W2 and the design rotative speed Nw, 
ulate by iteration with Vt3m obtained from ^ 


calcu 


''t3M “ ‘J 3M - ''Z3M ®3M 


and T3 - T2 obtained from 


(B12) 


(B13) 


^3 “ ^2 ■ ^^3 M*'^T3M ~ / ^2M^T2M^ 


the pressure coefficient corresponding to the input 9 from 


_ S^ad(^3 - T2) 


U 


3T 


(B14) 


(B15) 


with the following conditions; 

(1) W3 = Wg 

(2) Cp and y are functions of t3 obtained from subroutine CPF 

(^) 83M = 43Md a$3M<p, where A83M<p is obtained from an option 

to alter rotor deviation ahgle for off-design 9 values. 

frv»< calculations within CSPSI are repeated for every input 9 value 

for each stage of the compressor. ^ k y 


Subroutine CSPSD 

This subroutine is called by MAIN when the user has specified the option 
to alter the pressure rise coefficient 9 calculated at design speed for 
off-design speeds. Subroutine CSPSD calculates a change in the presLrr 
rise coefficient a*N for an off-design rotative speed N. 
calculation procedure is as follows: 

(1) Assume design values calculated in CSPREF for bomh. boma 
9d» and B3 |^jj. '*2Md» 

(2; For an input off-design rotative speed N. calculate V22M from 
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(B16) 


and the flow corresponding to the input N from 





(B17) 


(3) For this value of Wj; and the off-design rotative speed N, with 
Vi3(vi obtained from 



tan 63^ 


(B18) 


and T3 - t2 obtained from 



1— { U w 

Cp N^J^3M''T3M 


^2m'^T2M^ 



(B19) 


calculate by iteration the change in the pressure rise coefficient a+w 
corresponding to the off-design speed N from 








- V, 


tB20) 


with the following conditions: 

(1) W 3 = Wo 

(2) Cp and y are functions of t3 obtained from subroutine CPF. 

(3) 63IV1 03M(j + 63(vi|\j, where aB3mn is obtained from an option to 

alter the rotor deviation angle for off-design N values. 

These calculations are repeated within CSFSD for every input N value 
tor each stage of the compressor. CSP3D changes the pressure coefficient 
4' by an amount a4»|^ for an off-design part speed N. The overall 
etfect of a typical stage characteristic t(<p) is depicted in figure lb. 


Subroutine CSPAN 


CSPAN, which is called from MAIN, checks the value of input CB2M, CB2MR 
and CB3MR for each compressor stage. If any of this input is not equal to ’ 
^ero, a blade reset has been specified and CSPAN proceeds to alter the stage 
design flow coefficient 9(j and the pressure coefficient uint A 
new stage characteristic 4.(9) for the blade reset is calculated as 
tol lows: 

(1) Update the rotor inlet and o.^tlet flow angle from 
®2M ®2Md ^ ^®2M (B21) 
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Fiyure 16. - E«ccU of program opUon to alter pressure coefficient 
for off-design part speed on a typical stage characteristic ip 


I I ) 

^2Md 


■(B22) 


^3M “ ^3Md ^®3M 


(B23) 


(2) Assume design values calculated in CSPREF for other parameters. 
Calculate V;j^2M 


72M 


^2M 

tan + tan e 


Determine the change i 
reset aYq from 

'^22M 

= -07— - fti 

* f I 


2M 

n the flow coefficient 


(B24) 

a<p corresponding to blade 
^0 


(B2b) 



=1 


and the flow from 
- P2^'''z2M 


(B2b) 


(3) For this value of and the blade reset ayq, calculate by 
iteration the change in the pressure rise coefficient 
corresponding to blade reset &y from 


^p^ad,d^^3 

AH' = T"' 

^3T 


■ ^2^ 


- *, 


(Bc7) 
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with the following conditions: 


: * 


1) W3 ~ W9 

1 ) tp aru^ Y are functions of t^ and obtained from subroutine CPF, 

(3/ 03|v| - SjMd ^ ^^^3 My» where is obtained from ati option to 

alter rotor deviation angle because of blade reset. 

and where Vyjiv] is obtained from 


"t3M " '^3M - ’'Z3H «3M 


(B28) 


and 


T3 - T2 Is obtained from 


^ 3-^2 


= i- [U 


3m''t3M 


" ^2m'^T2M 


] 


(B29) 


These calculations are repeated within CSPAN for each stage of the 
como-essor. For an example of how CSPAN alters the stage characteristic 
1>(<p) toi blade reset, consider the case where the vane (or stator) just 
upstream of a stage is reset or rotated) by an amount AYn. For this 

^® 2 M = Aro ^® 2 M '‘nd A63M are both equal to zero Fioure 4 

tharthf characteristic will be altered by CSPAN. ^ Note 

fin! !^!fJ • f pressure rise coefficient t and the range of the 

flow coefficient 9 are both altered by the upstream stator reset ayq. 


Subroutine CSOUPT 

This subroutine calculates and writes individual stage and cumulative 

Pf'ameters for various selected speeds and flow 
conditions. Output written by CSOUPT is in either all SI or all U.S. 
cusi-omary units as specified by the value of the input oarameter UNITS. A 
summarv of the coding within CSOUPT is as follows: 

( 1 ) Read input, which consists of speed fraction N/N,j, lowest flow - 
low increment aW, and highest flow Wu. Calculations are 

fraction N/N^ at flow increments aW 

Ti om Wl to Wh. 

1 ( 2 ) Calculate the meanline representative velocity diagram at the rotor 
inlet mean line radius for flow W and Ap, where Cn and v are 
functions of t obtained from CPF. ° y <ue 

»u- Calculate the stage flow coefficient s = V79 m/U'>t Nonnallv 
this calculated 9 will be within the range of <p for t^is stage ^ 

'■■'•oni the stage characteristics for 

this stage, the stage pressure rise coefficient and stage adiabatic 

t na7;, Ho-ever if the calculated , for this stage ?s no! etthin 

this f l§w W tnd CSOUPT is coded to stop calculations for 

Uiis flow W and to print a message with the calculated 9 and a 

s atement that this stage is in a stall or choke condition. 

( ) Calculate the stage temperature and pressure ratio and cumulative 
compressor adiabatic efficiency and pressure ratio. cumulative 

(b) Calculate the meanline representative velocity diagram at the rotor 
outlet meanhne radius for flow W and A3, where C^ and v are 
functions of t obtained from CPF. p y dre 
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(6) Using the preceding calculated meanline velocity diagrams calculate 
(1) rotor and stator incidence angles and (2) rotor diffusion factor. 

(7) Write the preceding calculated stage and compressor parameters for 

the selected flow W and speed fraction N/N^. i 

This procedure is repeated for each selected flow for the speed 
fractions selected from the input speeds read by subroutine CSINPT. After 
all calculations are performed for a set of data for one compressor, 
subroutine CSOUPT is coded to return to the main routine MAIN at the 
beginning of MAIN, and another set of data for another compressor will be 
processed if it is available. 


Option to Alter q> Value for Off-Design Speeds 

Within the main routine MAIN there is an user option to alter the flow 
coefficient <p values for off-design speeds because of real flow 
effects. This option is executed for each stage according to the user's 
specified value of the input parameter SPDPHI. Figure 17 depicts the 
general effect of this program option to alter flow coefficient for 
off-design part speed on a typical stage characteristic 4>(<p)* The 
general effect is that the range of the <p values is increased for 
off-design part speed. The expression within MAIN used to alter <p is 
obtained'-from the following: 

* (i - ^) 


where 

A<P,^ = - <P 


(B31) 


and 
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(B32) 

Substituting equations (B32) and (B31) inti equation (B30) yields the 
altered for off-design part speed N as 




Equation (B33) is coded into MAIN to alter 


9 


(B33) 


for off-design part speed N. 


Option to Alter Rotor Deviation Angle 


Within subroutines CSPSI, CSPSD, and CSPAN there is an optional 
calculation that, if desired, alters the rotor deviation angle 6 r for 
off-design values of the flow coefficient 9 , rotative speed N, and blade 
setting angle respectively. This option will be executed for each 
stage, and the option is selected by means of the input parameters DRDEVP, 
DRDEVN, and DRDEV 6 for and a 6 Ryq* respectively. Figure 18 

shows the various angles associated with a typical rotor, mean line blade 
element. At the rotor outlet the relative flow angle S 3 m is related 
to the deviation angle by 


I 

^R * ®3M ■ '^3M 


(B34) 


So for a fixed rotor exit blade metal angle < 3 ^ 


Rotation 



Figure 18. - Angles associated with a typical rotor neanline blade element. 


31 




and the option to alter rotor deviation angle can be expressed in terms of a 
change in rotor exit relative flow angle Ae3w. The expression within 
subroutines CSPSI, JSPSD, and CSPAN used to alter rotor deviation angle is a; 

follows: 




where the subscript i represents <p, N, and yq for subroutines 
CSPSI, CSPSD, and CSPAN, respectively. The design relative velocity ratio 




is calculated within subroutine CSPREF. 
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EXAMPLE PROGRAM OUTPUT LISTING 
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The data set from the section Example Input Data Set was used as input 
data for the example program output listing given in this appendix. 
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APPENDIX D 

PROGRAM SOURCE CODE LISTING 


C 

c 

c 

c 

c 

c 


Kxx STGSTK - A COMPUTER CODE FOR PREDICTING MULTISTAGE AXIAL FLOW 
COMPRESSOR PERFORMANCE USING A MEANLINE STAGE STACKING METHO^ 

BY R. J, STEINKE, NASA LEWIS RESEARCH CENTER 

INTERMEDIATE^OUTPUT^''''^ MAGOR SUBROUTINES AND WRITES 

COMMON /VECTOR/ CPC0C6), TITLEC18), RT2(12), RH2(12), RT3(12), 

PSIREF(12), ETAREFC12), PHIDES(12, 9,8), 

' ^^ADES( 12 , 9 ,8 ) , PHI(12,9,8), PSI(12,9,8), ETA(12,9 
PETAa2), NSTAGEC12), PCTSPD<9), 

XBET2M( 12 , 9 ) , BL EIdc l2 , 9 ) , TT(13), PT(13), PR(12), TR(12), PR0(12), 
XTR0(12), ETA0(12), BET3MR{12,9), V22h(12:9), V23MC12,9): AREA2(12) 
X, AREA3(12), RM2C12), RM3C12), UT2(12), UT3C12), UM2(12), UM3(12) 
5,'5r2MRtl2,9),DPSISC12,9),RSOLM(12), ^2M(12), CB2mi2): CB2MR( 

I^I^*CM(12), RDFM(12), SK2M(12), S1NCM(12), BET3M(12 

GFlREF(12),ETAINPa2) 
DB2M(12,9), DB2MR(12,9),DB3Ma2,9),DB3MR 
X(12,9), B2NB3R( 12 , 9 ) , V3DV2RC12) ,DB3 mRG(12) 

X,DB3MRN(12,9), DB3MRR( 12,8) vou\ i,: j , UBinKfat 12 J 

COMMON -'SCALER/ RU, PI, G,AJ, RAD, RG, DCP, GJ, G2J, RPMRAD, NSTA 
■X, NSPE, NPTS, PO, TO, DESRPM, DESFLO, UNITS 

^dimension'hstg?12)°''^'®''^'^^°^^^'^'‘°''^^'°''°^''®' DRDEVN, DRDEVP 

DATA NSTG/1,2,3,4,5,6,7,8,9,10,11,12/ 

RU=1545 . 44 
PI=3. 14159217 
G= 32.1740 
AJ=778.12 
RAD=57. 29578 
DO 7 11=1,12 
7 NSTAGE(II)= NSTG(II) 

10 CALL CSINPT 

XXX CALCULATE FIXED PARAMETERS 

PO = POX144.0 

DO 20 1 = 1, NSTA = - 

AREA2(I)= (RT2(I)+RH2(I) )X(RT2(I)-RH2(I))XPI/164.0 
AREA3(I)= (RT3(I)4RH3(I))X(RT3(I)-RH3(I))XPI/144.o 
SQRT(RT 2 (I)XX 2 - AREA2C I )x?2 , 0/PI ) 

SQRT(RT3(I)XX2 - AREA3 (I ) X72 . 0/PI) 

RM2(I)xdESRPMxRPMRAD 
. RM3(I)XDESRPMXRPMRAD 
BET2M(I,1) = BET2M(I,1)/RAD 
RK2M(I) = RK2MCI) + CB2MR(I) 

SK2M(I) = SK2H(I) + CB2M(I + 1) 

CB2M(I)= CB2M(I)/RAD 

CB2MR(I)= CB2MRC D/RAD ^ 

CB3MR(I)= CB3MR(I)/RAD 
20 CONTINUE 
CALL CSPREF 

IF(ETADES(1,1,1) .EQ, 0.0) CALL CSETA 
IF (PSIDES(1,1,1).EQ.0.0) CALL CSPSI 
IF (SPDPSI.EQ.1.0) CALL CSPSD 
CALL CSPAN 
DO 51 1=1, NSTA 

BET2M(I,1) = BET2M(I,1) H RAD 
51 BET3MR(I,1) = BET3MR(I,1) X RAD 
C XXX WRITE INTERMEDIATE OUTPUT 
WRITE (6,2120) 

^ ° ^ FLOCAL(I,l) = FLOCAL(I,l)X0. 453592 

vT-1 UCTA^ SET2M(I,l),BET3MR(I,l),RINCN(I),RDFin(I),SINCM(I), 

IF (UNITS. EQ. 1.0) FLOCAL(I,l) = FLOCAL ( I , 1 )/0 . 453592 
2051 FORMAT (IlOH STAGE PHIREF PSIREF ETAREF DPHIA 

X DPSIA FLOCAL BET2M BET3MR RINCM RDFM,lbH S 


RM2(I)= 

RM3(I)= 

UM2(I); 

UM3(I)= 
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0006800 

0006900 

0007000 

0007100 

0007200 

0007300 

0007400 

0007500 

3007600 

0007700 

0007800 

0007900 

0008000 

0008100 

0008200 

0008300 

0008400 

0008500 

0008600 

0008700 

0008800 

0008900 

0009000 

0009100 

0009200 

0009300 

0009400 

0009500 

0009600 

0009700 

0009800 

0009900 

0010000 

0010100 

0010200 

0010300 

0010400 

0010500 

0010600 




ORIGINAL PAG'£ IS 
OF POOR QUALITY 

^nn J ' if ; ^ ^ ° ° 0 • PI 0 • 2 ) ) 

uu Dd I-*X;N5TA 

BET2M(I,1) = BET2M(I,1) / RAD 

1)?5 = UnIpe;“^ (NSTAGE(I).I = 1.12), (PCTSPD(J). (DPSIS(I,J). 1 = 1,12- 

^”^18S"??f SPDaia5?3X)//?l3Fl'4?) PABLE//40X, 13H STAGE NUP1BER/- 

DO 60 I=1,NSTA 
DO 60 J=1,NSPE 
DO 60 K=1,NPTS 

c DESIOH SPEEDS 

DO 70 1=1,NSTA 
WRITE (6,21201 

1F(HSPE.LT.4) GO TO 70 
WRITE (6,2120) 

IF(NSPE.LT,7) GO TO 70 
WRITE (6,2120) 

70 CONTINUE 

'"’SJtAGE N0,as..3<5X, 

20.4)) 

CALL CSOUPT 
GO TO 10 

2120 FORMAT (IHO////) 

END 


PSl 


ETA 


)/(9Fl- 


OOOOlOO 
0000200 
0000300 
0000400 
0000500 
0000600 
0000700 
0000800 
0000900 
0001000 
0001100 
0001200 
0001300 
0001400 
0001500 
0001600 
0001 700 
0001800 
000 1 900 
0002000 
0002100 
0002200 
0002300 
C002400 
0002500 
0002600 
0002700 
0002800 
0002900 
0003000 


SUBROUTINE CSINPT 

SUBROUTINE CSINPT READS AND WRITES THE INPUT 
COMMON /VECTOR/ CPCO(6), TITLE(18) RT2(l?l 
XRH3(12), PHIREF(12), PSIREF(il J ETAREF ll ' 
XPSIDES(12,9,8), ETADES(12.9,8), PHI(12;9.I ' 
X»8), DPHIAC12 )a DPSIA(12)> DETA(12) 

X 8 Et 2 H< 12 , 5 ), BlEED( 12 ;;t;’TffnS!^^f(l|;“?JJf 3 ; 


DATA 

RH2(12), RT3(i2), 
PHIDES(12,9,8), 
PSI(12,9,8), ETA(12,9- 
PCTSPD(9), 

TR(12), PR0(12),- 


X(12?9)^^B2MB3R(12^9f ®®2M(12,9), DB2MR ( 12 , 9 ) , DB3M( 12^9 ) , DB3MR 
X,5b3MRN(12?9)?'dB3MRP^ V3DV2R( 12 ) , DB3MRG( 12 ) 

x?P?e/np^ 

) READ?5?iJ00?END=999)'(T^ ^RDEVP 

1000 

^°°°118A4?J//)”^''^''^°^’^°” STACKING PROGRAM ////20X, 

SPEEDN, CHARTS, PO, TO, 

SPEEDN, CHARTS, PO, TO, 


10 


READ (5,1010) STAGES, 
WRITE(6,2010) STAGES, 
RG = RU/WTMOLE 
DCP = RG/AJ 
GJ = G«AJ 
G2J = GJ»‘2.0 
RPMRAD = PI/360.0 


WTMOLE,DESRPM,DESFLO 

WTMOLE,DESRPM,DESFLO 
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U U U J 1 u u 
0003200 
0003500 
0003^00 
0003500 
0003600 
0003700 
0003800 
0003900 
0 0 0 0 0 0 
OOO^ilOO 
000<j200 
000<f300 
000<^^00 
000^500 
000^600 
000^*700 
000^)800 
000-^900 
0005000 
0005100 
0005200 
0005300 
0Q05900 
0005500 
0005600 
0005700 
0005300 
0005900 
0006000 
0006100 
0006200 
0006300 
0006^00 
0006500 
0006600 
0006700 
0006300 
0006900 
0007000 
0007100 
0007200 
0007300 
0007^100 
0007500 
0007600 
0007700 
0007300 
0007900 
0008000 
0OC8100 
0008200 
0008300 
0008400 
0008500 
0008600 
0008700 
0008800 
0008900 
0009000 
0009100 
0009200 
0009300 
0009400 
0009500 
0009600 
0009700 
0009800 
0009900 
0010000 
0010100 
001 0200 
001 0300 




ORIGINAL PAGE IS 
OF POOR QUALITV 


X'jj.u ruRiiHi ^orxu.w/ 

2010 FORMAT (80H STAGES SPEEDS POINTS PO IN 

ILE WT DES RPM DES FLOW//(SF10 . 3)////) 

READ (5,1010) SPDPSI, SPDPHI, DRDEVG, DRDEVN, DRDEVP, 
WRITE (6,2011) SPDPSI, SPDPHI , DRDEVG, DRDEVN, DRDEVP, 

2011 FORMAT (60M SPDPSI SPDPHI DRDEVG DRDEVN 
XUNITS//(6F10.1)////) 

READ (5,1020) ( CPCO( I ) , 1=1 , 6 ) 

WRITE(6,2020) ( CPCO( I ) , 1=1 , 6 ) 

1020 FORMAT (3E20.8) 

2020 FORMAT (072H CPCO(l) CPC0(2) CPC0(3) 

1PC0(5) CPC0(6)//(6E12.5)////) 

NSTA = STAGEN 


TO IN 

UNITS 

UNITS 

DRDEVP 


CPCO(4) 


MO 


NSPE= SPEEDN 
NPTS = CHARTS 
READ (5,1010) 
READ (5,1010) 
READ (5,1010) 
READ (5,1010) 
READ (5,1010) 
READ (5,1010) 
READ (5,1010) 
READ (5,1010) 
READ (5,1010) 
READ (5,1010) 
READ (5,1010) 
READ (5,1010) 
READ(5,1010) 
WRITE(6,2030) 


(RT2CI), 1=1, NSTA) 
(RH2(I), 1=1, NSTA) 
(RT3(I), 1=1, NSTA) 
(RH3(I), 1=1, NSTA) 
(BET2M(I,1), 1=1, NSTA) 
(CB2M(I), 1=1, NSTA) 

(CB2MR(I), 1=1, NSTA) 


1=1, NSTA) 
1=1, NSTA) 
1=1, NSTA) 


(CB3MR(I), 

(RK2M (I), 

(RSOLM(I), 

(SK2M(I), 1=1, NSTA) 

(PR(I), 1=1, NSTA) 

(ETAINP(I), 1=1, NSTA) 

(NSTAGE(I),RT2(I),RH2(I),RT3(I),RH3(I),BET2M(I,1), 
XCB2M(I),CB2MR(I),CB3NR(I),RK2M(I),RS0LM(I),SK2M(I), 1=1, NSTA) 

2030 FORMAT (IlOH STAGE RT2 RH2 RT3 RH3 

X BET2M CB2M CB2MR CB3MR RK2M RSOLM,10H 

XSK2M//(5X,I5,4F10.A,5F10.2,F10.4,F10.2)) 

WRITE (6,2120) 

WRITE (6,2031) (NSTAGE(I), PR(1), ETAINP(I), 1=1, NSTA) 

2031 FORMAT (30H STAGE PR ETAINP//(5X, 15, 2F10 . 4) ) 

READ (5,1010) (PCTSPD(J), J=1,NSPE) 

2120 FORMAT (IHO////) 

READ (5,1010) (ETARAT(J),J=1,NSPE) 

WRITE (6,2120) 

WRITE (6,2121) 

2121 FORMAT (20H 
DO 21 1=1, NSTA 

READ (5,1010) (BLEEDd, J), J = 1,NSPE) 

21 CONTINUE 

WRITE (6,2120) 

WRITE (6,2041) 

1),J=1,NSPE) 


(PCTSPD(J),ETARAT(J), J=1,NSPE) 
PCTSPD ETARAT//(2F10 .4) ) 


(NSTAGE(I),I = 1,12), (PCTSPD(J), (BLEEDd, J), 1 = 1,12 


2041 FORMAT (20X,27H BLEED( STAGE, PCT SPD) TABLE//40X, 13H STAGE NUMBER/ 
18H PCT SPD,12(I5,3X)//d3F8.3)) 

DO 30 1=1, NSTA 

READ (5,1010) (PHIDES(I,1,K),K=1,NPTS) 

READ (5,1010) (PSlDESd,l.K),K = l,NPTS) 

READ (5,1010) (ETADESd,l,K),K = l,NPTS) 

30 CONTINUE 

DO 50 1=1, NSTA 

WRITE (6,2120) 

WRITE(6,2050) NST AGEd ) , (PCTSPDC J ) , J=1 , 3 ) , 

1,J,K),ETADES(I, J,K), J=1,3),K=1,HPTS) 

IF(NSPE.LT.4) GO TO 50 
WRITE (6,2120) 

WRITE(6,2050) NSTAGEd ), (PCTSPDC J ), J=4 , 6 ) , 

1,J,K) rcTADESd, J,K), J=4,6),K = 1,NPTS) 

IFCNSPE.lt. 7) GO TO 50 
WRITE (6,2120) 

WRITE(6,2050) NST AGEd ), (PCTSPDC J ), J = 7 , 9) , 

1, J,K) ,ETADESd, J,K),J = 7,9),K = 1,NPTS) 

50 CONTINUE 

2050 FORMAT (20X,41H INPUT DESIGN CHARACTERISTICS FOR STAGE— , I3//3(5X, 
1F10.3,10H PCT SPD ,5X)//3(30H PHIDES PSIDES ETADES)/(9F1 
20.4)) 

C CHANGE METRIC INPUT INTO ENGLISH UNITS 
IF (UNITS. NE. 1.0) GO TO 53 


((PHIDESd, J,K),PSIDESd 
((PHlDESd, J,K),PSIDESd 
((PHIDESd, J,K),PSIDESd 
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0010400 
0 0 1 0500 
0010600 

0 0 10700 
0010800 
0010900 
0011000 
0011100 
00 1 1200 
0011300 
0011400 

001 1500 
001 1600 
0011700 


0013500 

0013600 

0013700 

0013800 

0013900 

0014000 

0014100 

0014200 

0014300 

0014400 

0014500 

0014600 

0014700 

0014800 

0014900 

0015000 

0015100 

0015200 

0015300 

0015400 

0015500 

C015600 

0015700 

0015800 

0015900 

0016000 

0016100 

0016200 

0016300 

0016400 

0016500 

0016600 

0016700 

0016800 

0016900 

0017000 

0017100 

0017200 

0017300 

0017400 

0017500 

0017600 

0017700 

0017800 

0017900 

0018000 

0018100 

0018200 

0018300 

0018400 

0018500 

0018600 

0018700 

0018800 

0018900 












52 
51 

53 
999 


PO = PO/0.689<i76 
TO = T0K9. 0/5.0 
DESFLO = DESFLO/0. 953592 
DO 51 I=1,NSTA 
RT2(I) = RT2m/2.59 
RH2(I) = RH2(I)/2.59 
RT3(I) = RT3(I)/2.59 
RH3(I) = RH3(I)/2.59 


lyjllGINAL PAGE IS 
OF POOR QUALITY 


DO 52 J 
BLEEDCI, J) 
CONTINUE 
RETURN 
STOP 
END 


1 ,NSPE 
- BLEED(I,J)/0. 455592 


CALCULATES PARAMETERS AT DESIGN SPEED AND FLOW 


SUBROUTINE CSPREF 
: KKX SUBROUTINE CSPREF 
: CONDITIONS 

COMMON /VECTOR/ CPC0(6), TlTLE(lB), RT2(12)> RH2(12), RT3(12), 
XRH3n2), PHlREFa2), PS1REF(12), ETAREF(12), PHIDESC 12, 9i8 ) , 
XPSI.DEStl2,9,8), ETADES(12,9,8), PHI(12,9,8), PS1(12,9,8), ETA(12,9- 
X,8), DPHIA(12), DPSIA(12), DITA(12), NSTAGE(12), PCTSPD(9), 

TTC13), PT(13), PR(12), TR(12), PR0(12),- 
XTR0(12), ETA0(12), BET3MR(12,9), VZ2M(12,9), V23M(12,9), AREA2(12)- 
X, AREA3(12), RM2(12), RM3(12), UT2(12), UT3(12), UM2(12), UM3(12) - 
X,BET2MRa2,9),DPSIS(12,9),RS0LM(12), RK2M(12), CB2M(12), CB2MR( - 

X12), CB3MR(12), R1NCMC12), RDFM(12), SK2M(12), SlNCM(12), BET3M(12- 
X,9), PHIFIX(12), DPH1F(12),CPREF(12), GFlREFC 12 ) , ETAINP( 12) 
X,FL0CAL(12, 9) , ETARAT(9), DB2M(12,9), DB2MR(12 , 9 ) , DB3M( 12 , 9) ,DB3MR- 
X(12,9), B2MB3R(12,9), V3DV2R(12),DB3MRG(12)- 

X,DB3MRN(12,9), DB3MRP(12,8) 

COMMON /SCALER/ RU, PI, G, AJ, RAD, RG, DCP, GJ, G2J, RPMRAD, NSTA- 
X, NSPE, NPTS, PO, TO, DESRPM, DESFLO, UNITS 

X, CP, GAMMA, GM1,GF1,GF2,GF3,SPDPSI,SPDPHI,DRDEVG, DRDEVN, DRDEVP 
J *■ I 
1 = 1 

TT(I)= TO 
PT(I)= PO 

RHOT= PT(I)/(TT(I)xRG) 

TS= TT(I) 

RHOS= RHOT 

UT2(I)= RT2(I)xDESRPMxRPMRAD 
UT3(I)= RT3(I))iDESRPM>(RPMRAD 
CALCULATIONS AT ROTOR INLET 
VZ2M(I,J) = DESFL0/(RH0SXAREA2(I) ) 

V= VZ2M(I,J)/C03(BET2M(I,J)) 

CP= CPF(TS) 

RHOS= RHOTX(1.0-VKV/(G2JXCPKTT(I)))XXGF1 
TS= TT(I)H(RHOS/RHOT)xxGM1 
WCAL= RH0S>«AREA2(I)»«V22M(I, J) 

IF(ABS(WCAL - DESFLO) .GT. 0.01) GO TO 30 
CPREF(I) = CP 
GFIREF(I) = GFl 
PHIREF(I)= VZ2M(I,J)/UT2(I) 

DO 40 K=2,NPTS 

IF(PNIREF(I)-PHIDES(I, J,K)) 50,60,40 

CONTINUE 
K= NPTS 

PSIREF(I)= (PS IDES (I, J,K)-PSIDES(I, J,K-1))X(PHIREF(I)-PHIDES(I, J, - 
lK-1))/ (PHIDES(I,J,K)-PHIDES(I,J,K-D) + PSIDESC I , J , K*1 ) 
ETAREF(I)= (ETADESd, J,K)-ETADE5(I, J,K-1 ))x(PHIREF(I)-PHIDES( 1, J, - 
lK-1))/ (PHIDESCI, J,K)-PHIDES(I,J,K-1)) + ETADESd, J,K-1) 

GO TO 70 

PSlREFd)= PSIDESd,J,K) 

ETAREFd)= ETADESd, J,K) 

CONTINUE 

IF (P5IREFd).EG.O.O) GO TO 71 

PRd)= d.O + PSIREFd)»UT3d)xUT3d)/(GJ»CP«TTd)))xxGF2 
CONTINUE 

lF(ETAREF(i).EQ.0.0) ETAREFd)= ETAINPCI) 


20 


30 


40 

50 


60 

70 


71 
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0 0 1 9000 
0019100 
00 19200 
0019300 
0019<^00 
0019500 
0019600 
0019700 
0019300 
0019900 
0020000 
0020100 
0020200 
0020300 
0020600 
0020500 
OO2O6O0 C 
0020700 
0020800 
0020900 
0021000 
0021100 
0021200 
0021300 
0021600 
0021500 
0021600 
0021700 
0021800 
0021900 
0022000 
0022100 
0022200 
0022300 
0022600 
0022500 
0022600 
0022700 
0022800 
0022900 
0023000 
0023100 
0023200 


0023300 

0023600 

0023500 

0023600 

0023700 

0023800 

0023900 

0026000 

0026100 

0026200 

0026300 

0026600 

0026500 

0026600 

0026700 

0026800 

0026900 

0025000 

0025100 

0025200 

0025300 

0025600 

0025500 

0025600 

0025700 

0025800 


OF pOOH 


:auTY 


20 


TR(I)= 1.0 + (PR(I)><><GF3-1.0)/ETAREF(I) 

IF^(PSIREF(I)!eQ*!o!o) PSIREF(I) = GJXCPx(TT(I + l)-TT(I))XETAR£Fa) - 
X/UT3a)>‘>'2 
PT(i + n= pimxpRd) 

Tr<0(I)= TT(I + 1)/T0 
PR0(I)= PT(I+1)/P0 

1= in 

IF(I .LE. NSTA) go to 
DO 80 1=1, NSTA 

VT2M= VZ 2 MCI, J)XTAN(BET2rUI,J)) 

VT2tiR= UM2(I) - VT2M 
BET2NR(I,J)= ATAH2(VT2Fm,VZ2f1(I, J)) 

TS= TT(I+1) 

RH0T= PT(1+1)/(TT(I+1)xRG) 

RH0S= RHOT 

KXX CALCULATIONS AT ftOTOR OUTLET 
81 VZ3M(I,J)= DESFL0/(RH0SXAREA3(D) 

VT3M= (CpK(TT(i + U-TT<I) )XGJ + UM2( I ) XVT2M)/UM3 ( I ) 

V= SQRT(VZ3M(I, J)xx2 + VT3MH)(2) 

CP= CPF(TS) 

RHOS = RHOTX(1.0-VXV/(G2JXCPxTT(I+1)))XXGF1 
TS = TT(I+1)X(RH0S/RH0T)xx6M1 
WCAL= RIIOSxAREA3(I)XVZ3M(I, J) 

IF(ABS(WCAL-DBSFL0) .GT. 0.01) GO TO 
BET3M(I,J) = ATAN2(VT3M,VZ3M(I,J)) 

SINCM(I) = BET3M(I,J)><RAD - SK2t1CI) 

VT3MR= UM3(I) - Vt3M 
BET3MR(I, J)= ATAN2(VT3Mft,V23M(I,J)) 

RINCM(I)= BET2MR(I, J)XRAD - RK2M(I) 

V2MR= VZ2MC I , J )/COSC BET2MRC I , J ) ) 

V3MR= VZ3mi, J)/COS(BET3MR(I, J)) 

V3DV2R(I) = V3MR/V2MR 
RDFM(I)= 1.0 - V3MR/V2MR 
XRM2( I ) )/RSOLM( D/V2MR 
D62M(I,J) = BET2M(I,J) 

DB2MR(I,J) = BET2MR(I,J) 

DB3M(I,J) = BET3M(I,J) 

DB3MR(I,J) = BET3MR(I,J) 

B2MB3R(I,J) = BET2KR(1,J) 

80 CONTINUE 
100 RETURN 
END 


81 


+ (RM3a)xVT3M - RM2a)XVT2M)/(RM3a) + - 


BET3MR(I,J) 


RH2C12), RT3C12), 
PHIDES(12,9,8), 
PS1(12,9,8), EtA(12,9- 


FUNCTION CPF(TS) ^ 

: SUBROUTINE CPF(TS) CALCULATES SPECIFIC HEAT FROM STATIC 

: TEMPERATURE USING FIFTH DEGREE POLYNOMIAL 

COMMON /VECTOR/ CPC0(6), TITLEdS), RT2d2), 

XRH3d2), PHIF;EFd2), PSIREF(12). ETAREFd2), 

XPSlDESd2.9,8), ETADESC12.9,8), PHId2,9,8), ^ 

X,8), DPHIA(12), DPSIAd2), DETA(12), NSTAGE(12), PCTSPD(9), 
XBET2M(12,9), BLEEDd2,9), TT(13), PT(13), PR<12), TR<12), PRO(12),- 
XTR0d2), ETA0d2), BET3MRd2,9), VZ2M(12,9), VZ3Md2,9), AREA2d2)- 
X, AREA3(12), RM2(12), RM3d2)* UT2d2), UT3d2), UM2d2), UM3d2) - 
X.BET2MRd2,9),DPSISd2,9),RS0LMd2), RK2M(12), CB2M(12), CB?MR( - 
X12), Cr3MR(12), RINCMdZ), RDFMd2), SK2Md2), SlNCMdZ), BET3Md2- 
' X.9), PlIFIXdZ), DPHIFd2),CPREF(12), GF1REF(12) ,ETA1NPC12) 

X,FLOCAL(12,9), ETARAT(9), DB2Md2,9), DB2MRd2 , 9) , DB3Md2 , ^ , DB3MR- 
X(12,9), B2MB3Rd2,9), V3DV2R(12) ,DB3MRGd2)- 

X,0B3MRN(12,9), DB3MRPd2,8) . ..cta 

COMMON /SCALER/ RU, PI, G, AJ, RAD, RG, DCP, 6J , G2J, RPMRAD, NSTA- 
X, NSPE, NPTS, PO, TO, DESRPM, . DESFLO, UNITS 

X, CP, GAMMA, GMl,GFl,6F2,GF3,SPDPSI,SPDPHI,DRDEVG, DRDEVN, DRDEVP 
CPF = CPCOd )+(CPCO(2)+(CPCO(3)+(CPCO(A)+(CPCO(5)+ CPC0<6)XTS)XTS)- 
X XTS;XT$)XTS 
GAMMA = CPF/(CPF - DCP) 

GMl = GAMMA - 1.0 
GFl = 1.0/GMl 

GF2 = GAMMA/GMl 
GF3 = 1.0/GF2 
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00?5900 

oo;' 600 o 


0010700 

ooioaoo 

0010900 
0011000 
0011100 
0011200 
0011300 
OOIHOO 
0011500 
0011600 
0011700 
0011800 
0011900 
0012000 
0012100 
0012200 
0012300 
0012400 
0012500 
0012600 
0012700 
0012800 
0012900 
0013000 
0015100 
0U13200 
0013300 
0013400 
0013500 
0013600 
0013700 
0013800 
0013900 
0014000 
0014100 
0014200 
0014300 
0014400 
00 14500 
0014600 
0014700 


0011800 
0011900 
0012000 
0012100 
0012200 
0012300 
0012400 
0012500 
0012600 
0012700 
0012800 
0012900 
0013000 
0013100 
0013200 
0013300 
0013400 
001 3500 
0013600 
0013700 
00 13800 
0013900 


RETURN 

END 


OF PO « 


page 'S 
quauty 


generates adiabatic efficiency versus flow 


SUBROUTINE CSETA 
: XMX SUBROUTINE CSETA 
: COEFICIENT 

COMMON /VECTOR/ CPC0(6), TITLEdS), RT2(12), RH2(12), RT3(12), 
XRH3(12), PHIREF(12), PSIREF(12), ETAREF(I2), PHIDES( 12, 9,8) , 
XPSIDES(12,9,8) , ETADESC 12,9,8) , PHI(12,9,8), PSI(12,9,8), ETA(12^9 
X,8), DPHIA(12), DPSIA(12), DETAC12), N3TAGE(12), PCTSPD(9), 
XBET2M(12,9), BLEED(12,9), TT(13), PT(13), PR(12), TR(12), PR0(12), 
XIR0<12), ETA0(12), BET3MR( 12 , 9 ) , VZ2M(12,9), VZ3M(12,9), AREA2<12) 
X, AREA3(12), RM2(12), RM3(12), UT2(12), UT3(12), UM2(12), UM3<12) 
X,BET 2 MR(l 2 . 9 ),DPSIS( 12 , 9 ),RS 0 Lri( 12 ), RK 2 M( 12 ), C 62 M( 12 ), CB 2 MR( 

X12), CB3MR(12), RINCM(12), RDFM(12), SK2M(12). SINCM(12), BET3M<12 
X,9), PHIFIX(12), DPHIF(12),CPREF(12), GFIREFC 12 ) , ETAINP ( 12 ) 
X,FL0CAL(12,9), ETARATC9), DB2N(12,9), DB2MRC 12 , 9 ) , DB3M( 12 , 9 ) , DB3MR 
X(12,9), B2MB3RC 12, 9 ) , V3DV2RC 12 ) , DB3MRGC 12 ) 

X,DB3MRN(12,9), DB3MRPa2,8) 

G2J, RPMRAD, NSTA 

X, NSPE, NPTS, PO, TO, DESRPM, DESFLO, UNITS 

X,CP, GAMMA, GM1,GF1,GF2,GF3,SPDRSI,SPDPH1,DRDEVG, DRDEVN, DRDEVP 

) HXX TWO PARABOLAS ARE USED FROM STALL AND CHOKE TO DESIGN CONDITIONS 
DO 10 I-l.NSTA 

PSMPRS= (PHIDESd,.),!) - PHIREFd))xx2 


PHIREFd))xx2 


PCMPRS = (PHIDESd, J, NPTS) 

AS= -0 . ixETAREFd )/PSMPRS 
AC= -0.2 X ETAREFd )/PCMPRS 
BS= -2.0xPHIREFd)XAS 
BC= -2.0 xpHIREF(1)xaC 
CS= ETAREFd) + ASXPHIREFCI )XX2 
CC= ETAREFd) + ACXPHIREFCI ) XX2 
DO 20 K=l,NPTS 

IF(PHIDESd, J,K) - PHIREF(D) 11,12,13 

11 ETADESd,J,K) = (ASXPHIDESCI, J,K) + BS ) ^PHIDESd , J , K) + CS 
GO TO 20 

12 ETADESd,J,K)= ETAREFd) 

GO TO 20 

13 ETADE3< I J ,K) = ( ACxPHIDESd , J , K ) + BC ) XPHIDES (I , J , K ) + CC 
20 CONTINUE 

10 CONTINUE 
RETURN 
END 


SUBROUTINE CSPSI 

: XXX SUBROUTINE CSPSI CALCULATES PRESSURE COEFICIENTS FOR INPUT FLOW 
^ COEFICIENTS 

COMMON /VECTOR/ CPC0C6), TIfLE(18), RT2(12), RH2d2), RT3(12), 
XRH3d2), PHIREFd2), PSIREF(12), ETAREF(12), PHIDESd2 , 9, 8) . 
XPSIDES(12,9,8), ETADESd2,9,8), PHId2.9,8), PSI(12,9,8), ETA(12,9 
X,8), DPHIAd2), DPSIA(12), 0ETAd2), N5TAGE(12), PCTSPD(9), 
XBE^2M(12,9), BLEEDd2,9), TT(13), PT(13), PR(12), TR(12), PR0d2), 
XTRdl2), ETA0d2), BET3MR(12,9), V22M(12,9), VZ3Md2,9), AREA2d2) 
X, AREA3(12), RM2d2), RM3(12), UT2d2), UT3(12), UM2(12), UM3(12) 
X,BET2MRd2,9),DfSISd2,9),RS0LM(12), RK2M(12), CB2Md2), CB2MRC 

X12), CB3MRd2), RINCM(12), RDFM(12). SKZMdZ), SINCN(12), BET3Md2 
X,9), PHlFIXd2), DPHIFd2),CPREF(12), GFlREFd2) , ETAINPd2) 
X,FL0CAL(12,9), ETARAT(9), DB2M(12,9), DB2MRd2 , 9 ) , DB3M(12 , 9 ) , DB3MR 
Xd2,9), B2MB3R(12,9), V3DV2Rd2 ) , DB3MRGd2 ) 

X.DB3MRN(12,9), DB3MRP(12,8) 

GJ' G2J, RPMRAD, NSTA 

X. NSFE, NPTS, PO, TO, DESRPM, DESFLO, UNITS 

X^CP,GAMMA,GM1 ,GF1,GF2,GF3,SPDPSI ,SPDPHI ,DRDEVG, DRDEVN, DRDEVP 

DO 100 K^l.NPTS 
I-l 
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TT(I) = TO 
PT(1) = PO 

10 V?3M(I,J) = VZ2M(I,J)»«PHIDIES(I,J,K)/PHIREF(I) 

VL’MR = VZ.3M(I,J)/COS(nD2MR(I,l)) 

RET3flR(I,J) : DB3f1R(I,l) 

VTZM = VZ3M(I,J)HTAN(DB2M(I,1)) 

VZS = VIZMXHZ + VZ3M(I,J)*(X2 
RHOT = ri ( n/<TT( nxRG) 

RMOS = RHOTXd.O - V2S/(GZJXCPREF( I)?<TT( I ) ) JXXGFIREFC I ) 
CE5FI.C = RH0SxAREA2Cn»<VZ3M(I,J) 

TS = TT<n 
ID = 0 

11 VT3M = UM3(I) - VZ3ri(I, J)XTAN(BET3MR(I, J)) 

CP = CPF(rS) 

DT = <UM3(I)XVT3M - UM2( I)HVT2M)/(GJXCP) 

TSA = (DT + TT(I))/TT(I) 

PTA3 = PT(1)X(1.0 + ETADES<i,J,K)x(TRA - 1.0 ))><kGF2 
TTA3 = DT + TT<I) 

RHOT = PTA3/(TTA3XRG) 

V3S = VT3MXX2 + V7.3M( I , J ) XX2 

RHOS = RIIOTXCI.O - V3S/(G2JXCPXTTA3) )XXGF 1 

TS : TTA3X(RH0S/RH0T)XXGM1 

WCAL = RHOSxaREA3(1)KV23M(I,J) 

IF (TRA.GE.l.U) GO TO 12 
DT = 0.0 
GO TO 13 

12 ID = ID + 1 

V23M(I,J) = DESFLC/(RH0SxAREA3(I)) — ~ 

V3MR = VZ3M(I,J)/C0S(BET3MR(I,J)) 

OPTION TO ALTER ROTOR DEVIATION ANGLE FOR OFF DESIGN FLOW 

COEFICIENT 

IF (DRDEVP.EQ.1.0) 

XDB3MRP(1,K) = -( 10 . 00/RAD)x(\'3MR/V2MR - V3DV2R(D) 
BET3MR(I,J) = DB3MR(I,1) + DB3MRP(I,K) 

IF (ABS(WCAL-DESFLC).GT.O.Ol) GO TO 11 

13 PSIDESCI, J,K) = GJxCPxDTXETADESd, J,K)/UT3(I)XX2 
DV3DV2 = V3DV2Rd) 

FV3DV2 = V3MR/V2MR 
1 = 1 + 1 

IF d.LE.NSTA) GO TO 10 
100 CONTINUE 
RETURN 
END 
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SUBROUTINE CSPSD 

KXX SUBROUTINE CSPSD ALTERS PRESSURE RISE COEFICIENTS FOR OFF DESIGN \ 

SPEEDS i 

COMMON /VECTOR/ CPCOC6), TITlEdS), RT2(12). RH2(12), RT3(12i. - ! 

XRH3d2), PH1REF(12), PSIREFd2), ETAREF(12), PHIDES(12 , 9 » 8 ) , - \ 

XrSI0CSd2.9.8) , ETADES(12,9,8), PHId2,9,8). PSld2,9,8), ETA(12,9- \ 

X.8). DPHIAd2), DPSIA(12), DETA(12), NSTAGEdZ), PCTSPD(9)» - 1 

XBET2M(12.9). BLEEDd2.9), TT(13). PT(13), PRd2), IRd2). PR0(12),- 
XTR0(12), ETA0d2), BET3MR(12.9), VZ2Md2,9), VZ3M(12,9), AREA2d2)- 
X. AREA3(12), RM2(12), RM3(12), UT2d2), UT3d2), UM2(12)^ UM3(12) - \ 

X,BET2MRd2,9),DPSISd2,9),RSOLMd2), RK2M(12), CB2Md2), CB2MRC - j 

X12). CB3MR(12), RINCM(12). RDFM(12), SK2Md2), SINCM(12), BET3M(12- i 

X.9), PHIFJXdZ), DPHIFdZ) ,CPREFd2) , GFlREFd?) ,ETAINPd2) 

X,FL0CAL(12.9), ETARAT(9), DB2M(12,9), DB2MR (1 2 , 9 ) , DB3M(1 2 , 9) . DB3MR- 
X(12,9), B2Mn3R(12.9), V3DV2R(12 1 , DB3NRGd2 )- < 

X,DB3MRNd2,9), DB3MRP(12.8) ' 

COMMON /SCALER/ RU, PI, G, AJ, RAD, RG. DCP, GJ, G2J, RPMRAD, NSTA- 
X, NSPE, NPTS. PO. TO. DESRPM, DESFLQ, UNITS - J 

X, CP. GAMMA. GMl,GFdGF2,GF3,SPDPSI,SPDPHl,DRDEV6» DRDEVN.- DRDEVP ] 

DO 100 J = 1.NSPE ■( 

1 = 1 

TTCn = TO i 

PTd) = PO j 

10 VZiMd.J) = VZZMd.nxpCTSPDCJ) \ 

V.''riR = VZ3Md.jr/COS(DB2MRd.l)) 

RFTlMRd.JT = DBTMRd.U 


0ftl6!NAL PMt IS 
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0020900 

0021000 

0021100 

0021200 

0021300 

002UiOO 

0021500 

0021600 

0021700 

0021SOO 

0021900 

0022000 

0022100 

0022200 

0022300 

0022^t00 

0022500 
0022600 
0022700 
0022S00 
0022900 
0023000 
0023100 
0023200 
0023300 
0 023*'t00 
0023500 
0023600 
0023700 
0023800 
0023900 
002mOOO 
002''il00 
002^*200 
0020300 
0 0 26 •'♦0 0 


0016800 

0016900 

*(HH5000 

0015100 

0015200 

0015300 

0015600 

0015500 

0015600 

0015700 

001580C 

0015900 

0C160O0 

0016100 

0016200 

0016300 

0016600 

C016500 

0016600 

0 (J 1 7 0 0 
0016800 
0016900 
0 0 1 ; 0 0 0 
C017100 
0017200 

001 ;300 
C017600 
P 0 1 7 5 0 0 
0017600 
0017/00 
C 0 1 7 0 0 

r m / " c 0 
0018000 


ID = 0 

VT2M= V22[1(I,1)«PCT5PD( J)X TAN< DB2M( 1 . 1) ) 

V2S = \nzm^2 + VZ3f1(I,J)Ki<2 
RHOT = PT(n/(TT(nHRG) 

RHQS = RHOT^Cl.O - V2S/ ( G2 Ji^CPREFC I) XTT ( I ) ) )^K6F1REF( I ) 

DESFIC = RH05kaREA2(I)«VZ3M(I, J) 

TS = TT(I) 

11 VT3P1 = UM3(I)^PCTSPD( J) - '7Z3MC I . J ) )<TAN( BET3MR( I, J ) ) 

CP = CPFCTS) 

DT = (UM3(I))<VT3M - UH2( I )KVT. 71)/ (GJ1<CP)^PCTSPD( J) 

TRA = (DT + TT(I))/TT(I) 

PTA3 = PT(I)K(1.0 + ETAREF(I)H'(TRA ~ 1.0))1<^GF2 
TTA3 = DT + TT(I) 

RHOT = PTA3/(TTA3^.RG) 

V35 = VT3M>^x2 + VZ3M(I,J)«M2 

RHOS - RHOT)<(1.0 - V35/ ( G2 J ^CP^TT A3 ) ) 

T3 = TTA3k(RHOS/RHOT)i<HGM 1 
WCAL = RllOS^^AREA3(I)KVZ3M(I, J) 

IF (I.NE.l) GO TO 12 
DVZ3M = VZ3M(I,J) 

12 CONTINUE 

ID = ID + 1 

VZ3M(I,J) = DESFLC/(RH0S«AREA3(in 
V3f'lR = VZ3M(I>J)/C05(BET3hR(I, J)) 

OPTION TO ALTER ROTOR DEVIATION ANGLE FOR OFF DESIGN SPEEDS 
IF (DRDEVN.EQ.1.0) 

XDD3MRN(I,J) = -(10 .00/RAD)«(V3MR/V2MR - V3DV2R(D) 

BET3nR(I.J) = DB3MR(I,1) + DB3MRN(I,J) 

IF (ABS(WCAL-DESFLC) .GT.0.01) GO TO 11 

DPSIS(I,J) = GJkCPkDTHETAREF(I)/(UT3(I)^<PCTSPD( J) )K)<2 - PSIREF(I) 
DPSIS(I,J) = DPSIS(I,J) - DPSIS(I,1) 

1= I + 1 

IF (I.LE.NSTA) GO TO 10 
100 CONTINUE 
RETURN 
END 


SUBROUTINE CSPAN 

SUBROUTINE CSPAN ALTERS FLOW AND PRESSURE RISE COEFICIENTS FOR 
BLADE RESET 

COMMON /VECTOR/ CPCOC6), TITLE(18), RT2(12), RH2(12), RT3(12), 
XRM3(12), PHIREF(12), PS1REF(12), ETAREFC12), PHIDES ( 12 , 9 .8 ) , 
XPSIDES(12,9,8) , ETADES( 12, 9,8) , PHI(12,9,8). PSI(12,9,8), ETA(12,9- 
X,8), DPHIA(12), DPSIA(12>, DETA(12>, NSTAGE(12), PCTSPDC9), 
XBET2M(12,9), BLEED(12,9), TT(13), PT(13), PR(12), TR(12), PR0(12),- 
XTR0C12), ETA0(12), BET3MR ( 12 , 9 ) , V22M(12,9), V23mi2,9), AREA2C12)- 
X, AREA3(12), RM2(12), RM3(12)/ UT2(12), UT3(12), UM2(12), UM3(12) - 
X,BET2MR(12,9) ,DP5IS(12,9),R50LM(12), RK2M(12), CB2M(12), CB2MR( - 

X12), CB3MR(12), RINCM(12), RDFM(12), SK2M(12)/ SINCM(12)» BET3M(12- 
X,9), PHIFIX(12), DrHIF(12),CrREF(12), GFIREF ( 12 ) , ET AINP ( 1 2 ) 
X,FLOCAL(12,9), ETARAT(9), DB2M(12,9), DB2MR(12,9),DB3M(12,9).DB3MR- 
X( 12,9) , B2MB3R( 12, 9) , V3DV2R(12),DB3MRG(12)- 

X,DB3MRN(12,9), DB3MRP(12,8) 

COMMON /SCALER/ RU, PI, G AJ, RAD, RG, DCP, G J , G2J, RPMRAD, NSTA- 
X. NSPE, NPTS, FO, TO, DESRPM, DESFLO. UNITS 

X, CP, GAMMA, GMl ,Gri,GF2,GF3,SPDPSI ,SPDPHI,DRDEVG, DRDEV^I, DRDEVP 
J = l 
1 = 1 

TT(I)= TO 

PT(I)= PO 

90 TS= TT(1) 

DPHIA(I) = 0.0 
DFSIA(I) = 0.0 
DETAd) = 0.0 

ir((CB2MCI) + CB2MR(I) + CB3MRC I ) ) . EQ . 0 . 0 0 ) GO TO 93 
BET2M(1,J)= DB2M(I,J) + CD2M( I) 

BET2MR(I.J)= DB2MR(I.J) + CB2MR(I) 

BETTMRCI.J) = DB3MR(I,J> f CB3MR(I) 

VZrmi,J)= UM2(I)/(TAN(BET2M(I,J))+ T AN ( BET2MR ( I , J ) ) ) 

V2MR VZPM<i,J)/C05(BET2MRCT\J)) 
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0000100 

0000200 

0000300 

0000400 

0000500 

0000600 

0000700 

0000800 

0000900 

0001000 

0001100 

Q0C120O 

0001300 

0001400 

0001500 

0001600 

0001700 

0001800 

0001900 

0002000 

0002100 

0002200 

0002300 

0002400 

OOO^SOC 
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0002700 
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0002900 
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0003400 

0003500 
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‘ if -Jk-rt 7l 

C HHX CHANGE IN FLOW COEFICIENT FOR RESET 

DPH1A(I)= V22M(1, J)/UT2(1) - PHIREF(I) 

VZ3M(1,J)= V22mi,J) 

VT2M= V22M(I, J)?<TAN(DET2M(I,J)) 

V/2S = VT2NXH2 + V23ri( I , J ) KH2 
RMOT = PT(I)/<TT(UKRG) 

RHOS = RHOTxd.O - V2S/ ( G2J“CPREF( I )XTT (I ) ) )»«HGF1REF( I ) 
DE3FLC = RHOS>'AREA2(nKVZ3M(I,J) 

FLOCAL(I,J) = DESFLC 

92 VT3M=: UM3(I) - VZ3ri(I,J) xTAN( BET3MR(I , J ) ) 

CP = CPF(TS) 

DT= (UM3(I) X VT3M - UM2( I )XVT2M)/(GJKCP) 

TRA = (DT + TT<I))/TT(n 

PTA3 = PT(I)X(1.0 + ETAREF(I)X(TRA -1.0))XXGF2 

TTA3 = DT + TT(I) 

RHOT= PTA3 /(TTA3 KRG) 

V3S = VT3MXX2 + VZ3f1(i,J)XK2 

RHOS = RHOTXCI.O - V3S/(G2JXCPXTTA3 ))KXGF1 

TS = TTA3 X(RH0S/RH0T)XXGM1 

WCAL = RH0SXAREA3(1)xvZ3M(I, J) 

VZ3M(I,J) = DESFLC/(RH0SXAREA3(D) 

V3MR = VZ3M(I,J)/COS(BET3MR(I,J)) 

C XXX OPTION TO ALTER ROTOR DEVIATION ANGLE FOR BLADE RESET 
IF (DRDEVG.EQ.1.0) 

XDB3MRG(I) = -(10. 00/RAD)x(V3MR/V2MR - V3DV2R(i)) 
BET3MR(I,J) = DB3MR(I,J) + CB3MR(i) + DB3MRG(1) 

IF (ABS(WCAL-DESFLC) .GT.0.01) GO TO 92 
C XXX CHANGE IN PRESSURE RISE COEFICIENT FOR RESET 

DPSIA(I) = GJXCPXDT/(UT3(I)XUT3(I))XETAREF(I) - PSIREF(I) 

93 I = I+l 

IF (I.LE.NSTA) GO TO 90 
100 RETURN 
END 


SUBROUTINE CSOUPT 

C XXX SUBROUTINE CSOUPT CALCULATES AND WRITES STAGE AND COMPRESSOR 
C PERFORMANCE PARAMETERS 

COMMON /VECTOR/ CPC0(6), TITLEdS), RT2d2), RH2d2), RT3d2), 
XRH3d2), PKIREFd2), PSIREF(12). ETAREF(12), PHIDES(12 , 9 , 8) , 
XPSIDESd2,9,3), ETADESd2,9,8) , PHI(12,9,8), PSI(12,9,8), ETAd2,9 
X,8), DPHIA(12), DPSIA(12), DETA(12), NSTAGE(12), PCTSPD(9), 
XBET2M(12,9), DLEEDd2,9), TT(13), PT(13), PR(12), TRd2), PR0(12), 
XTR0(12), ETA0d2), BET3MR(12,9), VZ2M(12,9), VZ3M(12,9), AREA2(12) 
X, AREA3(12), RM2d2), RM3(12>, UT2d2), UT3(12), UM2(12), UM3d2) 
X,BET2MRd2,9),DPSISd2,9),RS0LMd2), RK2M(12), CB2M(12), CB2MR( 
X12), CB3MR(12), RINCM(12), RDFM(12), SK2Md2), SINCMdZ), BET3Md2 
X,9), PHIFIX(12), DPHIFd2),CPREF(12), GFlREFd2) , ETAINPd2) 
X,FL0CALd2,9), ETARAT<9), DB2Md2,9), DB2MRd2 , 9 ) , DB3M(12 , 9 ) , DB3MR 
Xd2,9), B2MB3Rd2,9), V3DV2R<12) , DB3MRG(12) 

X,DB3MRN(12,9), DB3MRP(12^8) 

COMMON /SCALER/ RU, PI, G, AJ, RAD, RG, DCP, GJ, G2J, RPMRAD, NSTA 
X, NSPE, NPTS, PO, TO, DESRPM, DESFLO, UNITS 

X, CP, GAMMA, GMl ,GF1 , GF2 , GF3 ,SPDPSI , SPDPHI , DRDEVG, DRDEVN, DRDEVP 

WRITE(6,2C70) 


2070 FORMAT (1H1////20X,27H XX COMPUTED OUTPUT DATA Mx////) 
C XXX READ SPEED AND FLOWS 

80 READ (5.1010) SPEEDF, FLOWIN,DFLOW, FLOWFl 
IF (UNITS. NE. 1.0) GO TO 81 

FLOWIN = FLOWIN/0. 453592 
DFIOW = DFLOW/0. 453592 
FLOWFl = FLOWFI/0, 453592 

81 CONTINUE 

IF (UNITS. EO. 1,0) FLOWIN = FLOWlNxQ . 453592 
WRITE(6,2080) SPEEDF. FLOWIN 
IF (UNITS. EQ. 1 . 0) FLOWIN = FLOWIN/O .453592 
2080 FORMAT (IHO//// 


1 17H PERCENT SPEED = , F8 . 3 , 1 OX, 1 5H INLET FLOW = ,F8.3//080 

IH STAGE PHI PSI TAU ETA PR 

2C-ETA C-PR.40H STG-FLOW R-INCM R-DFH S-INCM/) 

DO 82 J=1 ,NSPE 
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I 


I) 


k- 




0003700 IF (SPEEDF.EQ.PCTSPD(J)) JS=J 

0003800 82 CONTINUE 

U003900 C XHX CALCULATE THE OUTPUT 

OOOAOQO 1=1 

OOOAIOO TT(I)= TO 

000A20O PT<I)= PO 

OOOA300 WTFLOW = FLOWIN 

OOOAAOO 91 RHOT= PT ( I )/ ( TT ( I ) *<RG) 

000A500 TS= TT(I) 

O00A60O RHOS= RHOT 

0009700 C 

0009800 UTFLOW = WTFLOW - FLOWINXBLEED( I , JS) 

0009900 U2 = UT2(I)XSPEEDF 

0005000 U5 = UT3(I)«SPEEDF 

0005100 UMM2 = UM2 ( I ) XSPEEDF 

0005200 UMM3 = UM3( I ) ><SPEEDF 

0005300 100 V2 = WTFL0W/(RH0S»<AREA2(I)) 

0005900 V= VZ/C05(BET2M(I,in 

0005500 CP= CPF(TS) 

0005600 RHOS= RHOTX( 1 , 0-VXV/( G2J«CPXTT( 1 ) ) )x»Grl 

0005700 IF ((VXV) .GT. (G2J5«CPXTT(I))) GO TO 113 

0005800 TS= TT(I)H(RHOS/RHOT))««Gf11 

0005900 WCAL = RH0SXVZXAREA2( I) 

0006000 IF(ABS(WCAL-WTFLOW).GT.0.01) GO TO 100 

0006100 PHIC = VZ/U2 

0006200 IF (PHIFIX(I).NE.O.O) PHIC = PHIC + DPHIFfTl 

0006300 IF(PHIC.LT.PHI(I,JS,D) GO TO 110 

0006900 IF(PHIC.GT.PHI(I,JS,NPTS)) GO TO 120 

0006500 DO 130 K=2,NPTS 

0006600 IF(PHIC-PHI(I, JS,K)) 190,150,130 

0006700 130 CONTINUE 

0006800 K= NPTS 

0006900 190 PSIC = (PSI(I, JS,K)-PSI(I, JS,K-1))x(PHIC-PHI(I, JS,K-1) )/ 

0007000 1 (PHKI, JS,K)-PHI(I,JS,K-D) + PSI ( I , JS , K-U 

0007100 ETAC=(ETA(I, JS,K)-ETA(I, JS,K-1))X(PHIC-PHI(I, JS,K-in/ 

0007200 1 (PHKI, JS,K)-PHI(I, JS,K-D) + ET A ( I , JS , K-l ) 

0007300 GO TO 160 

0007900 150 PSIC= PSI(I,JS,K) 

0007500 ETAC= ETA(I,JS,K) 

0007600 160 CONTINUE 

0007700 PR(I) = (1.0 + PSIC»U3«U3/ (GJ«CP)«TT( I > ) )KXGF2 

0007800 TAU = PSIC/ETAC 

0007900 TR(I)= 1.0 + ( PR ( I ) «KGF3-1 . 0 )/ETAC 

OOOSOOO TT(I+1)= TT(I)hTR(I) 

0008100 PT(I+1)= PT(I) XPR(I) 

0003200 TRO(I)= TT(I+l)/TO 

0008300 FRO(I)= PT(I+l)/PO 

0008900 IF(I.EQ.l) GF3S = GF3 

0008500 GF30= (GF3 + GF3S)/2.0 

0008600 ETAO(I)= (PRO( I )iO«GF30 - 1.0)/(TRO(I) - 1.0) 

0003700 VT2M = VZ X TAN ( BET2M( I , 1 ) ) 

0003800 VT2MR = UMM2 - VT2M 

0003900 BET2MR(I, JS)= ATAN2(VT2MR, VZ) 

0009000 RINCM(l): BET2MR(1,JS) « RAD - RKZMd) 

0009100 V2MR= VZ/COS( BET2HR( I , JS) ) 

0009200 RHOT= PT ( K 1 )/ ( T T ( I + l) >(RG) 

0009300 TS= TT(I+1) 

0009900 RHOS= RHOT 

0009500 161 VZ3M(I,JS)=WTFL0W/(RH0SXAREA3(I)) 

0009600 VT3M= ( CPx ( TT ( 1 +1 )-TT ( I) ) xGJ + UMM2 XVT2M)/UMM3 

00097 00 VS= VZ3M(1,JS)hx 2 + VT3M>*x2 

0009800 CP= CPF(TS) 

0009900 RHOS= RHOT* ( 1 . 0-VS/ ( G2 J>*CP«TT ( I + l ) ) ) MMGFl 

0010000 IF((VS) .GT. (G2J*CP*TT(I+1))) GO TO 113 

001 0100 TS= TTd+l) * (RHOS/RHOT)**Gril 

0010200 WCALi RII0S*AREA.3d)*VZ3Md, JS) 

0010300 IF(ABS(WCAL-WIFLOW).GT.0.01) GO TO 161 

001 0900 BET3nd,JS) = ATAN2(VT3M,VZ3Md,JS)) 

0010500 SlNCMd) = BET3Md,JS)><RAD - SKZMd ) 

0010600 VT3MR = UMM3 - VT3M 

0010700 BET3MRd,JS)= ATAN2( VT3MR . VZ3Md , JS ) ) 

0010800 V3MR = VZ3Md , JS )/CCS ( BET3MRd , JS) ) 

001 0900 RDFMd): 1.0 - V3MR/V2MR ♦ (RM3d)*VT3M - RM2d )*VT2M)/(RM3d) + 
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0011000 

0011100 

0011200 

0011300 

ooinoo 

0011500 

0011600 

0011700 

0011300 

0011900 

0012000 

0012100 

0012200 

0012300 

0012400 

0012500 

0012600 

0012700 

0012S0G 

0012900 

0013000 

0013100 

0013200 

0013300 

0013400 


XRM2(I) )/RS0LM(I)/V2MR 
IF (UNITS. EQ. 1.0) WTFLOW = WTFLOW*0 . 453592 
C WRITE THE OUTPUT 

IF (UNITS. EQ. 1 .0) WTFLOW = WTFLOW/0 . 453592 
2090 FORMAT ( 5X, 15 , 8F10 . 4, FIO . 2 , FIO . 4, FI 0 . 2) 

1= I + l 

IF(I.LE.NSTA) CO TO 91 
GO TO 111 

110 WRITE(6,2100) I,PHIC 
GO TO 111 

120 WR1TE(6, 2110)1, PHIC 
GO TO 113 


2100 FORMAT (lOH FOR STAGE, 13, 18H , 
2110 FORMAT (lOH FOR STAGE, 13, 18H , 

111 FIOUIN = FLOWIN + DFLOW 

IF (FLOWIN. LE.FLOWFI) GO TO 81 
113 IF (JS.lt. NSPE) GO TO 80 
DO 112 I=1,NSTA 
DO 112 J=1,NSPE 

112 DPSISd.J) = O.O 
RETURN 

1010 FORMAT (8F10.0) 

END 


COMPUTED PHI IS,F8.4,06H STALL) 
COMPUTED PHI IS,F8.4,06H CHOKE) 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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